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In this work we aim at understanding the effect of n-type and p-type substitutional doping in the 
case of matrix-embedded and freestanding Si nanocrystals. By means of ab initio calculations we 
identify the preferential positioning of the dopants and its effect on the structural properties with 
respect to the undoped case. Subsequently, we consider the case of phosphorus and boron co-doped 
nanocrystals showing that, against the single-doping situation, the energetics strongly favors the 
binding of the impurities at the NC surface. Indeed we demonstrate that the polar B-P bond forms 
a stable permanent electric dipole that radially points inwards the nanocrystal. Such noteworthy 
characteristic and its physical consequences are discussed alongside new potential applications. 


I. INTRODUCTION 

The use of silicon nanocrystals (Si-NCs), embedded 
in wide-gap materials, or in colloidal solution, has 
found large employment in applications such as opto¬ 
electronics and photonics?® photovoltaics?®^ an d bio- 
imaging???^ Several advancements in the fabrication 
techniques permit nowadays to control with good pre¬ 
cision the size, shape, and density of the Si-NCs. Col¬ 
loidal Si-NCs allows, in principle, the fabrication of 
large-area optoelectronic devices by vacuum-free print¬ 
able processes. Usually, in order to avoid NCs agglom¬ 
eration one functionalizes the nanocrystals surface with 
organic molecules, but this can be detrimental to the 
electron transport properties? 4 Thus, the production of 
colloidal Si-NCs without surface functionalization pro¬ 
cesses is highly desirable. This has been recently ob¬ 
tained by synthesizing all-inorganic Si-NCs using simul¬ 
taneous doping with n- and p-type impurities? 815 4 ' Re¬ 
garding matrix-embedded Si-NCs, they offers advantages 
in terms of stability, low-cost manufacturability and the 
development of CMOS compatible devices. However, the 
high resistivity of the embedding matrix reduce the pos¬ 
sibility to employ such composite material in applications 
that require a high carrier mobility. One way to enhance 
the current of Si-NCs-based devices is to increase the 
density of the NCs in the samples, thus increasing the 
silicon content and reducing the tunneling distance be¬ 
tween the Si-NCs. The downside of this approach lies 
in the difficulty of controlling the density of the Si-NCs 
distinctly from other parameters (e.g. size, shape, and 
so on), and in the limited efficiency of quantum confine¬ 
ment while approaching the percolation threshold, thus 
affecting the optical properties? 8 In alternative, the pos¬ 
sibility of introducing dopants has been explored in re¬ 
cent years in order to improve transport and other funda¬ 
mental characteristics?^^ Clearly, in the case of Si02- 
embedded Si-NCs one may focus on doping the sole Si02 
matrix, thus improving the conductivity of the sample 
while not altering the optical absorption of the Si-NCs, 


the latter tunable for example by changing the NCs size 
distribution. 22 

Differently than bulk materials, the doping of nanostruc¬ 
tures involves additional complications arising from the 
presence of interfaces, and of materials with different 
affinities in hosting the dopant species?^ In practice, ex¬ 
perimental samples of embedded Si-NCs present inhomo¬ 
geneous impurity concentrations, that depend on the fab¬ 
rication technique!® In other words, since dopants tend 
to diffuse during the fabrication process, one cannot sim¬ 
ply limit the doping to a selected region of the material. 
Therefore, in order to correlate the effect of the dopants 
to the observed quantities, it is important to understand 
the preferential positioning of the impurities in the nanos¬ 
tructures. In this context, the presence of defects in 
unpassivated samples will likely favor the presence of 
dopants in the proximity of interfaces (i.e. interstitial 
doping) J® In practice, at low impurity concentrations 
the impurity atoms will act as passivants, correspond¬ 
ingly improving the optical yield. Conversely, in the case 
of high concentrations or of well passivated (or undefec¬ 
tive) samples, the impurities substitute to the host ma¬ 
terial atoms, by diffusing towards the most energetically 
favorable sites!® 

In the following we investigate the effect of n-type and 
p-type substitutional doping in the case of Si-NCs em¬ 
bedded in Si02 and of freestanding, colloidal Si-NCs. By 
means of ah initio total-energy calculations we identify 
the preferential location of dopants in the nanostructure, 
and its effect on the structural properties with respect to 
the undoped case. Finally, the case of co-doped NCs will 
be investigated and discussed, in particular for colloidal 
NCs, alongside its potential applications. 


II. STRUCTURES AND METHOD 

All the structural and electronic properties have been 
obtained by full ab-initio calculations in the frame¬ 
work of density functional theory (DFT) using the 
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FIG. 1: Ball-and-stick representation of a SiC> 2 -embedded Si-NC with an impurity atom placed inside the NC (left panel) and 
at two different interfacial sites, corresponding to Si 3+ (central panel) and Si 1+ (right panel) in the undoped structure. Red 
and cyan spheres represent O and Si atoms, respectively. Si-NC atoms are connected by thick sticks, while the yellow region 
marks the position of the dopant. 


SIESTA package! 25 * 26 ! Calculations have been performed 
using norm-conserving pseudopotentials within the local- 
density approximation (LDA). An energy cutoff of 250 Ry 
on the electron density and no additional external pres¬ 
sure or stress have been applied. Consistently with con- 
vergency trial runs, all the calculations have been per¬ 
formed with a Gamma point k-sampling and standard 
double-^ basis set for all the atoms. Atomic positions 
and cell parameters have been left totally free to move. 
The SiC> 2 -embedded NCs have been generated starting 
from a 3x3x3 betacristobalite-SiCE matrix, by removing 
all the oxygen atoms inside a sphere whose radius deter¬ 
mines the NC size. The so-obtained Si-NC, embedded in 
the SiC> 2 , presents perfectly coordinated atoms and the 
same Si-Si distance of betacristobalite, corresponding to 
about 3.2 A. After the relaxation, the Si/SiC >2 interface 
forms strongly stressed bonds, while far from the inter¬ 
face the bulk atomic densities are recovered.^ 

We have introduced the dopant impurity in a substitu¬ 
tional Si-site located in the NC center, at the Si/SiC >2 
interface, or in the SiC >2 far from the NC. Since the inter¬ 
facial Si atoms can neighbour with one-to-three oxygen 
atoms (i.e. Si 1+ , Si 2+ , Si 3+ ), for each structure we have 
considered all the possible oxidation type for the interfa¬ 
cial dopingPS Each structure has been fully relaxed after 
doping in order to include the effect of the structural 
rearrangements around the dopant atom. Examples of 
so-obtained doped structures are shown in Figure [1] 

All the calculations have been performed for the Si 32 , the 
Si 47 , and the Si 35 NCs embedded in SiC> 2 . However, since 
we have obtained equivalent results for the three different 
cases, for the sake of simplicity we will present only the 
results of the latter system in the following. 

The freestanding colloidal Si-NCs have been obtained 


from bulk silicon by a spherical cutoff, and by passi¬ 
vating all the dangling bonds with OH groups (oxidized 
Si-NCs). The choice of OH as Si-NCs terminatio n has 
several reasons. First of all, we have already shown 27 2Sl 
that, beside strain effects, there is a strong similarity be¬ 
tween the calculated electronic and optical properties for 
the silica embedded Si-NCs and for the Si-NCs simply 
passivated by the OH groups. Thus these properties are 
almost completely determined by the barrier provided by 
the first shell of oxygen atoms. More over i t is worth to 
note that, from the present and othePMl calculations, 
the energetics of the doping sites seems crucially affected 
only by the termination type while the presence of an 
embedding matrix should play a minor role. Like for the 
embedded Si-NCs, in our calculations each freestanding 
structure has been fully relaxed after doping. We have 
used a 3nm-sized periodic cubic cell for the freestand¬ 
ing NCs smaller than the Sis 7 (diameter ~ 1.36nm), and 
5 nm for the others. The error on the total dipole due to 
the interaction with the NC replica has been estimated 
by applying, for the largest NC, the Martyna-Tuckerman 
correction to the scf potential in order to mimick the 
case of a truly insulated NCf 3 -^ The so-obtained total 
dipole was about 6% smaller than the uncorrected one; 
this value should indicate the maximum error on the cal¬ 
culated dipole moments. 

In the following we will consider the Sis 7 -(OH) 76 as ref¬ 
erence system for the study of co-doping. 
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III. SINGLE-DOPED SILICON 
NANOCRYSTALS 

We report in Figure [2] the total energy comparison for 
the Si 35 NC (diameter ~ 1 nm) embedded in Si02, single- 
doped with a group-V (N or P) or a group-III (A1 or B) 
atom placed at different substitutional sites. The choice 
of taking into account N impurities is related to the fact 
that N presents very different electronegativity, oxida¬ 
tion states, and electron affinity, with respect to P. Con¬ 
versely, As atoms should behave similarly to P, given an 
almost exact matching of the above characteristics. 

In Figure [2] the minimum of each curve identifies the 
energetically favored site of the corresponding dopant. 
The data clearly indicate that for n-type doping (P 
and N) the impurity tends to settle in the NC core, 
while for p-type doping (A1 and B) the interfacial sites 
are favored. These results are on line with recent the¬ 
oretical calc ulation s made on oxidized freestanding Si 
nanoparticlesP^ 3 1 Moreover they agree excellently with 
the three-dimensional atomic probe tomography experi¬ 
mental data of Duguay and coworkeri 21 322 that demon¬ 
strated that n-type dopants are efficiently introduced in 
the ’’bulk” of the Si-NCs, whereas B atoms are preferen¬ 
tially located at their periphery, at the Si/SiC >2 interface. 
Furthermore, we note from Figure [2] that in all the con¬ 
sidered cases the doping of the SiC >2 region is unlikely 
to occur. In particular, for N and P doping is required 
a very high formation energy to move the dopant atom 
from the NC core to the silica (3.9 eV and 7.9 eV re¬ 
spectively), in agreement with the observation that for 
Si nanoclusters embedded in a SiC >2 matrix, the matrix 
provides a strong barrier to P diffusion inducing P seg¬ 
regation in the Si rich regionPS Instead, still consistently 
with experimental result^!, the diffusion of B and A1 
toward SiC >2 results significantly easier, with barriers of 
just 1.1 eV and 1.6 eV, respectively. 

Actually, the presence of defects in real samples should 
foster interstitial doping in the silica and at the Si/Si 02 
interface. Such condition should be superposed to 
the above one in which substitutional doping acts on 
the bas is of the lowest formation energy in defect-free 
samples ! 20 ! 24 ^ 

It is worthwhile to note that a similar picture has been 
observed for doped Si nanowires (NWs). Fukata et 
alP*l, in a joint experimental-theoretical effort, have been 
shown that, whereas the impurity B atoms are more sta¬ 
ble at the Si-SiC >2 interface in a split interstitial configu¬ 
ration at a Si site, P impurities preferentially piles up in 
the crystalline Si region. 

The correspondence of the dopant positioning in different 
Si nanostructures, independently on their size or shape, 
indicates that the energetics is mainly governed by the 
atomic binding energy of the dopant, while other param¬ 
eters (e.g. strain) play a minor role. 

To better understand the effect of the different dopants 
on the NCs structural properties we have reported in Fig¬ 
ure [3] the length of the bonds formed by the dopant with 


the four neighbouring atoms of the substituted Si. We ob¬ 
serve that P and A1 maintain a sort of tetra-coordination 
in all the cases, with all the bonds showing approximately 
constant characteristic lengths with Si and O. Instead, N 
tends to form in all the case a strong anti-bonding char¬ 
acteristics with one of the neigbours, that is repulsed 
at large distance. The same occurs for B, but only when 
placed at the interface with two oxygens (interface2), cor¬ 
responding to the lowest-energy configuration. 

To further investigate the local configuration at the im¬ 
purity site, we provide in Table [T] the charge population 
analysis of all the structures, obtained by the Mulliken 
method- 36 and the electronegativities of the considered 
atomic species. 

Electronegativity plays a crucial role in determining 
the polarity of the bonds, i.e. the amount of charge lo¬ 
calized onto a certain atom. It is interesting to note that 
the most stable sites correspond for the n-type dopants 
to the highest values in terms of population, while for 
the p-type dopants to the lowest ones. Due to the much 
higher electronegativity of oxygen with respect to silicon, 
the charge population of the dopant decrease when mov¬ 
ing from the NC core to the silica. Such expectation is 
fulfilled in all cases except when the P and B dopants 
are localized in the silica. In our opinion this should be 
related to the proximity of the four neighboring oxygen 
atoms (see Fig. [3] top-right panel), leading to an overes- 
timation of the charge population. 



core interfacel interface2 silica 

FIG. 2: Total energy of the SEs NC embedded in the Si02 
structure and doped with N, P, Al, or B, substitutional of a 
Si atom in the NC center (core), at the interface bonded to 
one or two oxygens (interfacel, interface2), or in the SiC >2 far 
from the NC (silica). For a better comprehensibility each line 
has been shifted so that the structure with the dopant in the 
NC center has zero energy. 
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FIG. 3: Length of the bonds formed by the dopant with its 
neighbours, for different dopant site in the Si 35 NC embedded 
in SiC >2 (please refer to Fig. [ 2 ] for nomenclature). Circles 
highlight the most stable site for each dopant. 


Mulliken population 


Si 35 /Si0 2 

core 

interfacel interface2 

silica 

N 

5.38 

5.12 

5.00 

4.83 

P 

5.13 

4.99 

4.86 

5.08 

Al 

3.73 

3.27 

2.84 

2.57 

B 

3.31 

3.29 

3.20 

3.39 


Electronegativity 


B 

C 

N 

O 

2.04 

2.55 

3.04 

3.44 

Al 

Si 

P 

S 

1.61 

1.90 

2.19 

2.58 


TABLE I: Top: Mulliken population of the dopant atom for 
the Si 35 NC embedded in Si02 and doped with either N, P, 
Al, or B, substitutional of a Si atom in the NC center (core), 
at the interface bonded to one or two oxygens (interface 1, 
interface2), or in the Si0 2 far from the NC (silica). Bold text 
highlights the most stable site for each dopant (see Fig. [2|. 
Bottom: Pauling electronegativities of the two upper lines of 
the periodic table, from group III to VI. 


IV. CO-DOPED SILICON NANOCRYSTALS 


Coming to the case of multiple dopant atoms, in par¬ 
ticular to that of compensated co-doping, one should cor¬ 
respondingly increase the size of the hosting Si-NCs, in 
order to avoid a forced aggregation of the impurities. As 
already stated above, the energetics of the doping sites is 
mainly affected by the termination type at the interface, 
while the presence of an embedding matrix plays only 
a minor role. Therefore we have employed freestanding 
oxidized NCs as an approximation of the fully embed¬ 
ded ones, significantly reducing the computational effort 
and making possibile to consider larger NCs, that are 
more suitable in the presence of multiple dopant atoms. 


Moreover this choice allows us to take contact with the 
experiments performed on all-inorganic colloidal Si-NCs , 
co-doped with B and P, and terminated by OH ! 13 I 16 I 1,j 

The case of co-doped nanostructures presents interest¬ 
ing aspects related to the possibility of the dopant atoms 
to interact and bind together. In particular, in the case of 
opposite dopant types (i.e. n- and p-type), the emergence 
of new chemical pathways that can alter the picture ob¬ 
tained in the sing le-d o ped case has been proved for both 
Si NCs and NWs.™®l 

To shed light on these aspect we have considered, as start¬ 
ing point, a freestanding oxidized NC, Sis 7 (OH) 7 6 , co¬ 
doped with P and B atoms located at their energetically 
favored sites, as determined above for the single-doped 
case: P at the NC center and B at the NC surface. Fur¬ 
ther, to investigate the effect of the distance and of the 
interaction between the dopants we have considered addi¬ 
tional nanostructures in which the P atom progressively 
approaches the B atom (see Fig. [4j left panel). We have 
numbered the position of the P atom from 5 (center) to 
1 (inner surface). From the analysis of the total energy, 
reported in Figure [5] it is clear that while in the single- 
doped NC the P atom stabilize at the NC center, in the 
co-doped case the formation of a strong P-B bond (po¬ 
sition 1) is clearly energetically favored, with a gain of 
about 1.4 eV. 

In alternative, the P-B bond can be formed at the NC 
center, fulfilling the energetics of the P dopant in the 
single-doped case while forcing the B atom in the NC 
core (see Fig.[4j right panel). Such configuration, named 
lc, results more energetically favored than the 5, 4, 3, 
and 2 ones, confirming the idea that the formation of the 
P-B bond rules over the other parameters (see Fig. |5|. 
Hovewer, the P-B pair at the interface (position 1) re¬ 
veals more stable than that at the NC center (position 
lc), with a difference in the formation energies of about 
0.3 eV . 

Due to the polar nature of the P-B bond, the above con¬ 
dition implies the formation of a static electric dipole ra¬ 
dially directed and preferentially located at the NC sur¬ 
face, pointing toward the NC center. The high stability 
of the dipole direction is evidenced by considering a con¬ 
figuration in which the co-dopants are switched, named 
li (see Fig. [4j center panel). Quite interestingly, such 
configuration produce a total energy higher than all the 
other systems, and about 1.45 eV higher than the most 
stable case. 

We have computed the induced dipole as the difference 
between the total electric dipole of co-doped and undoped 
NCs. For the most stable configuration we have calcu¬ 
lated an induced dipole moment /i ~ ID, corresponding 
to two opposite charges of about 0.1 a.u. placed at a dis¬ 
tance of 2.32 A (the P-B bond length). 

Despite the modest magnitude of the single dipole, the 
simultaneous presence of many inward-pointing dipoles 
shall be sustained by the strong stability of the dipole 
direction. In the case of the Si 8 7 -(OH) 76 doped with six 
P-B pairs disposed along the three directions (two pairs 
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FIG. 4: Freestanding OH-terminated Sis 7 NC (left panel) doped with B at the surface and P placed at one of the sites from 
the center (5) to the nearest neighbour of B (1); (central panel) doped with P placed at the surface and B in the inner nearby 
(li); (right panel) doped with P at the NC center and B in the outer nearby (lc). Cyan, red, white, green, and yellow spheres 
represent Si, O, H, B, and P atoms, respectively. 


for each axis), we have computed a difference in the for¬ 
mation energy of 28 meV when switching one of the P-B 
pair atom positions (i.e. inverting the dipole direction), 
meaning that up to six inward-pointing dipoles are ener¬ 
getically permitted at environmental temperature at this 
NC size. Besides, we expect that a much greater num¬ 
ber of dipoles should be supported in a large NC, due to 
larger dipole-dipole internal separations (i.e. lower repul¬ 
sive energies). The latter consideration is in agreement 
with Sugimoto and coworkers, reporting six dopant pairs 
in colloidal NCs of diameter 1.5 nm, and 10-1000 dopant 
pairs in NCs of diameter 2-10 nm, respectively.^ 

In a series of paperjh ^ 16 * 17 !, Sugimoto and cowork¬ 
ers suggested that in the strong doping regime the co¬ 
doped NCs should be depicted by an internal shell posi¬ 
tively charged surrounded by an external shell negatively 
charged, implying substantial repulsive forces between 
neighboring NCs that enable the stability of the colloidal 
Si-NCs in polar liquids without the formation of precip- 



FIG. 5: Total energy of the co-doped freestanding Sis 7 -(OH )76 
as a function of the dopants position (see Fig. |4j). Inset: total 
dipole moment of the undoped NCs, terminated by H or OH 
groups, as a function of the NC size. 


itates. 

However, unless some charge is injected during the NCs 
preparation, we can restrict to the case of neutral NCs. 
In such case, since no electric field is emitted by two (ho¬ 
mogeneous) concentric spherical shells of opposite charge, 
we would expect vanishing repulsive forces for large NCs 
in the strong co-doping regime (i.e. homogeneous charge 
limit). 

We worth to stress that, differently than simple hy¬ 
drogenated NCs, oxygenated NCs show a non-negligible 
“macroscopic” dipole due to a strong charge trapping at 
the surface oxygens.®^ To make this clear, we have per¬ 
formed a total-dipole calculation on a set of fully-relaxed 
NCs of different size, H- and OH-terminated, reported in 
the inset of Figure [5] 

Therefore, while our calculations support the shell model 
proposed by Sugimoto et al., we have a different belief 
about the mechanism of stability. In our opinion, the sur¬ 
face dipoles form to the detriment of the “macroscopic” 
dipole of the NC, naturally existing in the undoped case. 
As a consequence, while the undoped NCs aggregate by 
the “macroscopic” dipoles interaction, the neutralization 
of such dipoles by co-doping produce unpolarized NCs 
that cannot aggregate. Eventually, the screening pro¬ 
vided by the polar liquid, and the thermal energy, may 
help in contrasting the aggregation induced by a residual 
polarization. 

This view is supported by our case of NC with six co¬ 
dopants, presenting a total dipole about 7% smaller with 
respect to the undoped case. Clearly, for larger NCs con¬ 
taining hundreds of co-dopants, we expect a much larger 
dipole-cancellation effect. Beside we also note that the 
undoped Sig 7 NC shows the largest dipole of the set (see 
Fig. [5] inset), thus difficult to cancel out. 

In the case of co-doped embedded NCs, the surface 
dipoles could in principle be aligned by an external field 
during the annealing, and would remain so afterwards 
due to the cage effect of the embedding matrix. If so, the 
possibility of forming a silicon-based ferroelectric semi- 
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conductor may constitute a convenient alternative to con¬ 
ventional materials made by toxic or rare elements.^ 


V. CONCLUSIONS 

By means of theoretical calculations we have substitu- 
tionally positioned a dopant atom inside or outside Si- 
NCs embedded in a Si02 matrix, revealing that whereas 
p-type (A1 and B) dopants preferentially site at the NC 
interface, the n-type ones (N and P) are energetically 
favored in the Si-NC core, in agreement with recent ex¬ 
perimental observations! 21 * 24 -^ The energetics of the sys¬ 
tems as a function of the dopant position shows that Si 02 
forms a very large diffusion barrier for P and N, and a 
much reduced one for B and Al. The analysis of the bond 
lengths and charge population around the dopant atom 
reveals a correlation between the electronegativity of the 
impurity and its preferential location, a strong indica¬ 
tion that the presented results should be generalizable to 
other Si nanostructures like oxygenated NWs or Si/SiC >2 
superlattices. 

Following the results for the single-doped Si-NCs, in the 
case of n-type and p-type compensated co-doping we have 
investigated the most stable configuration for P+B im¬ 


purities. We have revealed that the formation of a P- 
B bond is energetically favored, with B positioned at 
the NC surface and P at the inner nearby. As a con¬ 
sequence, NC co-doping involves the spontaneous for¬ 
mation of static dipoles at the NC surface, very stable 
and radially directed inwards the NC. In our view, the 
role of these dipoles is that of disrupting the “macro¬ 
scopic” dipole of the NC which is associated to aggrega¬ 
tive forces in colloidal suspensions, and can provide an 
alternative explaination about the possibility of produc¬ 
ing stable colloidal Si-NCs without any surface func¬ 
tionalization process.^ The presented outcomes can be 
very compelling for a large set of applications, from 
photovoltaics j 5 6I 10 I 43 I to luminescence^, bioimagingj 12 ! 
and nanomechanicsamong others. 
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